Time-optimal synthesis of unitary transformations 
in coupled fast and slow qubit system 
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In this paper, we study time-optimal control problems related to system of two coupled qubits 
where the time scales involved in performing unitary transformations on each qubit are significantly 
different. In particular, we address the case where unitary transformations produced by evolutions 
of the coupling take much longer time as compared to the time required to produce unitary trans- 
formations on the first qubit but much shorter time as compared to the time to produce unitary 
transformations on the second qubit. We present a canonical decomposition of SU(4) in terms of the 
subgroup SU(2) x SU(2) x U(l), which is natural in understanding the time-optimal control problem 
of such a coupled qubit system with significantly different time scales. A typical setting involves 
dynamics of a coupled electron-nuclear spin system in pulsed electron paramagnetic resonance ex- 
periments at high fields. Using the proposed canonical decomposition, we give time-optimal control 
algorithms to synthesize various unitary transformations of interest in coherent spectroscopy and 
quantum information processing. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 



The synthesis of unitary transformations using time- 
efficient control algorithms is a well studied problem 
in quantum information processing and coherent spec- 
troscopy. Time-efficient control algorithms can reduce 
decoherence effects in experimental realizations, and the 
study of such control algorithms is related to the com- 
plexity of quantum algorithms (see, e.g., [3, [!,[!])• Signif- 
icant literature in this subject treat the case where uni- 
tary transformations on single qubits take negligible time 
compared to transformations interacting between differ- 
ent qubits. This particular assumption is very realistic 
for nuclear spins in nuclear magnetic resonance (NMR) 
spectroscopy. Under this assumption Ref. [j (see also 

B S 0, SI, M EL HI 0, HI EE El EE El) presents 
time-optimal control algorithms to synthesize arbitrary 
unitary transformations on a system of two qubits. Fur- 
ther p rogr ess in the case of multiple qubits is reported in 



ing the significantly different time scales immanent in the 
system. Employing this canonical decomposition, we de- 
rive time-optimal control algorithms to synthesize various 
unitary transformations. Our methods are applicable to 
coupled electron-nuclear spin systems occurring in pulsed 
electron paramagnetic resonance (EPR) experiments at 
high fields, where the Rabi frequency of the electron is 
much larger than the hyperfine coupling which is fur- 
ther much larger than the Rabi frequency of the nucleus. 
In the context of quantum computing similar electron- 
nuclear spin systems appear in the Refs . l2rl l27l . [H, l29l . 

IMlILIIIliBSlHlIi 

The main results of this paper are as follows. Let 5^ 
and I v represent spin operators for the fast (electron spin) 
and slow (nuclear spin) qubit, respectively. Any unitary 
transformation G G SU(4) on the coupled spin system 
can be decomposed as 



ther p rogr ess in the case of multi 

n Ha ee ei m m m m m i 



G = K\ expihS 13 !, + t 2 S a I x )K 2 , 



(1) 



In this work, we consider a coupled qubit system where 
local unitary transformations on the first qubit take sig- 
nificantly less time than local transformations on the sec- 
ond one. In addition, we assume that the coupling evo- 
lution is much slower than transformations on the first 
qubit but much faster than transformations on the sec- 
ond one. We present a canonical decomposition of SU(4) 
in terms of the subgroup SU(2) x SU(2) x U(l) reflect- 
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where S a I x and S°I X correspond to x-rotations of the 
slow qubit, conditioned, respectively, on the up or down 
state of the fast qubit. The elements K\ and K 2 are rota- 
tions synthesized by rapid manipulations of the fast qubit 
in conjunction with the evolution of the natural Hamil- 
tonian. The elements K\ and K 2 belong to the subgroup 
SU(2) x SU(2) x U(l), and in appropriately chosen ba- 
sis correspond to block-diagonal special unitary matrices 
with 2 x 2-dimensional blocks of unitary matrices. 

The minimum time to produce any unitary transfor- 
mation G is the smallest value of (\t\ \ + \t 2 \)/u>{., where u>l 
is the maximum achievable Rabi frequency of the nucleus 
and (ti, t 2 ) T is a pair satisfying Eq. (fTJ). Synthesizing K\ 
and K 2 takes negligible time on the time scale governed 
by 



2 



The paper is organized as follows. In Sec. HH we recall 
the physical details of our model system exemplified by a 
coupled electron-nuclear spin system. The Lie-algebraic 
structure of our model is described in Sec. IIII1 which 
is used to derive control algorithms (pulse sequences) 
for synthesizing arbitrary unitary transformations in our 
coupled spin system. In Sec. IIV1 we present examples. 
We prove the time-optimality of our control algorithms 
in Sec.fVl and some details of the proof are given in Ap- 
pendix [AJ 

Our work draws some results from the theory of Lie 
groups, which are explained as needed. We refer to [43|, 
|44| for general reference. To make the paper broadly 
accessible, we work with explicit matrix representations 
of Lie groups and Lie algebras. 



II. PHYSICAL MODEL 

As our model system, we consider two coupled qubits. 
We introduce the operators 5 M and l v which correspond 
to operators on the first and second qubit, respectively. 
In particular, these operators are defined by 5 M = (07, ® 
id 2 )/2 and I v = (id 2 ®cr v )/2 ( see [45| ) , where u x : — (} J), 
<jy := (? "o*)) and o z := (0-1) are the Pauli matrices 
and id 2 := (01) is the 2 x 2-dimensional identity matrix. 
In the remaining text, let /z, v £ {x, y, z} and 7 S {x, y}. 

In an experimental setting using an electron-nuclear 
spin system, the first qubit is represented by the electron 
spin (of spin 1/2). Similarly, the second qubit is rep- 
resented by the nuclear spin (of spin 1/2). We assume 
that in the presence of a static magnetic field pointing in 
the z-direction, the free evolution is governed in the lab 
frame by a Hamiltonian of the form 

i^ ab = u s S, + ujl, + J(2SJ Z ), (2) 

where ujs and cjj represents the natural precession fre- 
quency of, respectively, the first qubit and second qubit 
and J is the coupling strength. We assume that 

w s »w/» J. (3) 

This assumption is motivated by coupled electron-nuclear 
spin system occurring in EPR experiments at high fields 
(see, e.g., Sect. 3.5 of [46j]). The time scales in Eq. ([3]) 
insure that the hyperfine coupling Hamiltonian between 
the spins averages to the Ising Hamiltonian 2S Z I Z , as in 
Eq. @. This is the so-called high field limit. 

The first and second qubit are controlled by transverse 
oscillating fields, which result in the corresponding con- 
trol Hamiltonian given by H l g h + H\ &h , where 

H l s ah = 2w?(t) cos^fi + 4>s{t)]S x 
is the control Hamiltonian of the first qubit and 



is the control Hamiltonian of the second qubit. The 
amplitude, frequency, and phase of the control func- 
tion w.r.t. the first qubit are represented by ui^(t), wf , 
and 4>s = 4>s{t) respectively. Similarly, 0Jj.it), and 
fa = 4>i{t) represents the amplitude, frequency, and 
phase of the control function w.r.t. the second qubit. We 
use uj^. and u>f to denote the maximal possible values of 
u>*(t) and u>r(t). In our model system, we assume that 

W r 7 «J« CO? . (5) 

Therefore, we refer to the first qubit as the fast qubit and 
the second qubit as the slow qubit. 

We choose = u>s and u/ c = ojj — J. In a double 
rotating frame, rotating with the first and second qubit 
at frequency u>f and u>^, the transformations J7i a b(t) and 
U m t(t) describe, respectively, a unitary transformation in 
the lab frame and the double rotating frame. We have 

Uub(t) = exp(-itUcS z )exp(-itull z )U Tot (t). 

Using the rotating wave approximation, the Hamiltoni- 
ans i?o ab , H l g° , and H l f h transform, respectively, to 

H =JI Z + J(2S Z I Z ), (6) 
cos (f>s (t) + Sy sin (j>g (t)] , 

and 

Hi = loI (t) [I x cos cj)i (t) + I y sin fa (t)] . 

In absence of any irradiation on qubits, the system 
evolves under the free Hamiltonian —iHo- From the time 
scales in Eq. ([5]), we can synthesize any unitary transfor- 
mation of the form exp(— itS^) in arbitrarily small time 
as compared to the evolution under Hq or Hq + Hj. 
Let us define the operators, 

S * = (id 4 /2 + 

and 

where id<j is the d x c?-dimensional identity matrix and O2 
is the 2 x 2-dimensional zero matrix. Note that Hq = 
2JS@I Z , and the system is described by the Hamiltonian 

H + Hj = 2JS fl I z +w I r (t)(S a + SP)(I x cos 4>i+I y sin fa). 

Since J S> W^.(t), and S^Ij, does not commute with 
S^Iz, the above Hamiltonian gets in the first order ap- 
proximation truncated to 

H a (fa) = 2JSPl z +w I r (t)S a (I x cosfa + I y smfa). (7) 

Similarly, we can prepare an Hamiltonian 



Hf h = 2c^(t) cos[c4* + fa{t)\I x 



(4) 



H {fa) = 2JS a I z + wl{t)S^{I x cos fa + I y sin fa) (8) 
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a)i — J 



a(3 



lL>S + J 



cos ~ J 




FIG. 1: The eigenstates of the Hamiltonian _ffg ab are shown, 
where the transitions aa «-> a/3 and /3a «-> (3(3 correspond re- 
spectively to the orientation along and opposite to the static 
magnetic field. The first and second index refer to the orien- 
tation of the electron and nuclear spin, respectively. Refer to 
the text for details. 



by using H^(<j)i) — exp(iTrS x )H a ((j)j) exp(— inS x ). 

The Hamiltonians H a ((f>i) and H^((f>j), operate on the 
slow qubit and induce transitions aa <-> a/3 and f3a <-> (3(3 
of the nuclear spin as shown in Fig. Q] (cp. Table 6.1.1 of 
|4^|). The a and (3 states of the spins denote their ori- 
entation along and opposite to the static magnetic field, 
respectively. For the electron spin, the (3 state has lower 
energy than the a state as its gyromagnetic ratio is neg- 
ative. Similarly, for the nuclear spin, the a state has 
lower energy than the (3 state as its gyromagnetic ratio 
is positive (as for a proton). We remark that the energy 
eigenstates (3a, (3(3, aa, and a(3 correspond, respectively, 
to the basis states 00, 01, 10, and 11. In Fig. [TJ the first 
and second index in eigenstates refers to the orientation 
of the electron and nuclear spin, respectively. In absence 
of any irradiation on the two qubits, the system evolves 
under the Hamiltonian —iHq. In this section, we have 



shown how to synthesize generators of the form 
and -iH Q . 



-iS u 



ILL LIE-ALGEBRAIC STRUCTURE OF THE 
MODEL SYSTEM 



All transformations of our model system are contained 
in the Lie group G — SU(4), which is the set of 4x4- 
dimensional unitary transformations of determinant one. 
The operators —il^, ~iS v , and —i2I^S u , are infinitesi- 
mal generators of the Lie group G, and they generate the 
15-dimensional Lie algebra g = su(4) given by the (real) 
vector space of 4x4-dimensional (traceless) skew Hermi- 
tian matrices. We have shown how to synthesize gen- 
erators of the form —iS^, —iH a (4>j), and —iHo. These 
generators are sufficient to produce any unitary transfor- 
mation on the coupled qubit system, as described below. 



Lemma 1. The Lie algebra generated by the elements 
— iSp, —iH a {(f)[), and —iH , is equal to q — su(4). 

Therefore, a standard result on the controllability of 
(Thm. 7.1 of Ref. [4lJ) implies that the system is com- 
pletely controllable and any unitary transformation in 
G = SU(4), can be synthesized by alternate evolution 
under the above Hamiltonians. 

Lemma 2. The Lie algebra t, generated by the elements 
—iSfj, and —iHo consists of the elements —iS^, —i2S v L z , 
and — il z . 

The Lie algebra 6 represents a class of generators that 
take significantly less time to be synthesized, as they only 
involve controlled rotations of the fast qubit and evolu- 
tion of the free Hamiltonian —iHq (no controlled rota- 
tions of the slow qubit are involved). We can decompose 



(9) 



where the subspace p (of g) consists of the elements — ^ 
and —i2S t j j L 1 . The decomposition of Eq. ([9]) is a Cartan 
decomposition (see, e.g., [43j, p. 213) as 



[t,l] C 6, [t,p] C p, and [p,p] C f, 



(10) 



where [51,32] = .91.92 — .9231 is the commutator (5^ £ g). 

Let K = exp(t) denote the subgroup of G = SU(4) 
which is infinitcsimally generated by t. The elements 
of K can be synthesized only by the free evolution and 
employing controlled transformations on the fast qubit. 
Therefore, synthesizing transformations of K takes signif- 
icantly less time as compared to general unitary trans- 
formations not contained in K. In particular, controlled 
transformations on the slow qubit are necessary to syn- 
thesize general unitary transformations. The Lie group 
K = exp(6) is equal to S[U(2) x U(2)], which is sometimes 
referred as SU(2) x SU(2) x U(l). 

Consider a maximal Abelian subalgebra a contained 
in p. In our case, o is spanned by the operators —iS"l x 
and —iS a L x . Any element a S a can be represented as 
ai(—iS l3 I x )+a2(—iS a L x ), where at, 02 G M. As a matrix, 
a takes the form 



/0 01 0' 

a x 

a 2 

\0 a 2 0. 



We obtain the Lie group A = exp(a) corresponding to 
the Abelian algebra a. From a Cartan decomposition 
of a real semisimple Lie-algebra as satisfying Eqs. ©- 
(fT0|) . we obtain a dec omp osition of the compact Lie group 
G = KAK (see, e.g., [H, Chap. V, Thm. 6.7): 

Lemma 3. Any element G £ SU(4), can be written as 



G = K X eMti{^S p I x ) + t 2 {-iS a I x )]K 2 , 



(11) 



where ti,t% € R and K l7 K 2 E K. 

Remark 1. The computation of KAK decompositions 
was analyzed in Refs (ia. [49l. |50|. l5ll. |52| . In this work, we 
consider the Cartan decomposition, which corresponds 
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to the type AIII in the classification of possible Cartan 
decompositions (see, e.g., pp. 451-452 of Ref. [43jl. 

Transforming all elements G G G to SWAP ■ G ■ SWAP, 
where 



R 2 Ri (*■)_, Ro 



SWAP = exp(-wrS-I) 



/l 0\ 

10 

10 

\0 1/ 



S - {S x ,S y ,S z ) T , and I = {I x ,I y ,I z ) T , the KAK de- 
composition is given in explicit matrices by 



U x 2 
2 U 2 



exp 



Ui 2 

o 2 u 2 





/0 


a i 














a 2 






ai 





°i 











0/ 




C\ 


- 


-isi 








C-2 





-is 2 


— is\ 





Cl 








—is 2 





C 2 



t/ 3 o 2 
o 2 c/ 4 



c/ 3 o 2 
o 2 c/ 4 



where Sj = sin(aj/2) and Cj = cos(aj/2). In particular, 
the Lie group K is given in this basis by block-diagonal 
unitary transformations, where 2 is the 2 x 2-dimensional 
zero matrix and U\, U 2 (and C/3, U4) are 2 x 2-dimensional 
unitary matrices such that the product of their deter- 
minants is one. The considered KAK decomposition is 
equivalent to the cosine-sine decomposition [53|, HH, [5f| . 

Remark 2. In Ref. [4j, a different Cartan decompo- 
sition is considered. In that case, the subalgebra t is 
given by the elements — iS^ and — il v and corresponds 
to unitary transformations on single qubits of a coupled 
two-qubit system. Synthesizing unitary transformations 
on single qubits is assumed in Ref. [H to take significantly 
less time, as compared to unitary transformations which 
interact between different qubits. 

Since elements of K can be synthesized in negligible 
time, we obtain as the main result of this paper that the 
minimum time to synthesize any element G G SU(4) is 
the minimum value of (|ii | + |t 2 |)/cj^ such that [t\, t 2 ) T is 
a pair satisfying Eq. (fTTj) . We defer the proof of this fact 
to Sec. El Let us describe how to use the KAK decom- 
position of G, to synthesize an arbitrary transformation 
using only the generators —iS^, —iH a {<f>i), and —iHq. 

The Lie algebra 6 decomposes to ?i © pi, where ti is 
a subalgebra, composed of operators —iS^ and —i2S v I z , 
and pi is generated by —il z which commutes with all ele- 
ments of Ji . The Lie algebra ?i can be further subdivided 
by a Cartan decomposition ti = 6 2 ©p 2 . The subalgebra 
t 2 is generated by the operators —iS^, and the subspace 
p2 consists of the operators —i2S^L z . Therefore, similar 
as in Lemma 3, we obtain a decomposition of K: 

Lemma 4- Each element Kj G K can be decomposed 
as Kj = exp(-iT 2 j-iI z )L 2: j-iexp(-iT 2 j2S z I z )L 2 j = 

exp[-i(r 2j _i -T2j)I z ]L 2 j-iexp(-iT2jHo/J)L 2 j, (12) 

where Tj G M. and Lj G K 2 = exp(t 2 ). 
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FIG. 2: The figure shows a canonical pulse sequence for syn- 
thesizing unitary transformations in the coupled qubit system. 
Let R2 = R2 exp(iirS x ) and -Ro = exp(— ivoS z ) exp(— inS x ). 
Since 1/J ^ l/w r , the length of the time intervals tj/w r is 
larger as depicted. Refer to the text for details. 



Using an Euler angle decomposition (see, e.g., pp. 454- 
455 of Ref. [56[), the elements Lj G K 2 are given as 

Lj = cxp(-i9 jA S z ) expi-iOj.^S.j;) exp(~i6j, 3 S z ) 
= exp[-i(0 jtl + jt3 )S z ] exp[-iej t2 R(0 jt3 )}, (13) 

where R(9j, 3 ) = S x cos ^,3 — S y sin 6^,3. 

Similarly, any element A of the subgroup A can be 
written as A = exp[ii (-iS^Ix) + t 2 (-iS a L x )} = 



exp 



in 1 



e it 3 I, e -iUH Q /J exp 



-i^H a (0) 

in* 



e it 3 h e -it 1 H 13 (t 3 ) /w' r e -itiH / J e -tt 2 H a (Q) /wl 



(14) 



for £3 = 2Jt\/w I r mod 4-7T and £4 = J{t\ — t 2 )/w^ 
mod 27T > 0. This follows by substituting for expressions 
of H , H 01 ^), and ff^i) (see Eqs. ©-©). Com- 
bining Eqs. (fl"2")) - (fT4")) . a complete decomposition of an 
element G G SU(4), can be written as K\AK 2 = 



e -touS Be -iwi. Rie -inHo/J e X p 



xe -iUHo/J exp 



-i^H<*(r) 

in 1 



w r 

R 3 e- 1TiHo/J R 4 , 



where all the transformations Rj operate on the fast 
qubit. In particular, we have i?4 = exp[— i9^ 2 R{9^ 3 )\ 1 
R3 = exp[-i6^ 2 R(v 3 )], R 2 = exp[-i6 2t2 R(v 2 )], Ri = 
exp[-^i i2 i?(vi)], v 3 = 6» 3i3 + 6» 4i i + 6*4,3, v 2 = 9 2:3 + 9 3A + 
«3, «i = #1,3 + #2,1 + v 2 , v a = 6*1,1 + Vi, r = T 4 -r 3j and 
w = T\ — t 2 + t 3 — T4 — t 3 . The time to produce G is 
essentially {t\ +t 2 )/wj,. Note that exp(— iwl z ) = 



exp[-iwH /(2J)}e lnb *exp[-iwH /(2J)}. 



Transformations on the fast qubit such as exp(— ivoS z ) 
are significantly faster. Figure [2] shows the canonical 
pulse sequence realizing any unitary transformation as a 
sequence of rotations under —iHg, —iH^^j), and —iS^. 



5 



IV. EXAMPLES 



(a) 



Wc introduce the unitary transformations CNOT[l, 2], 
CNOT[2, 1], and SWAP which are given as follows 



/l 







n 












/l 




1 


















, and 


1 


















1 




1 




\0 


1 


0/ 


Vo 


1 





0/ 




\0 





Let c G {1,3,-1,-3}. The elements of SU(4) corre- 
sponding to the transformation CNOT[2, 1] are given 
by exp[c7r(— i2S x I z + iS x + iI z )/2), which is equal to 
exp(ic7r/4)CNOT[2,l]. For CNOT[l,2] and SWAP wc 
obtain the elements exp[c7r(— i2S z I x + iS z + iI x )/2] and 
exp[c7r(j2S , x / 2 ; + i2S y I y + i2S z I z ) /2], which are equal 
to exp(ic7r/4)CNOT[l,2] and exp(ic7r/4)SWAP, respec- 
tively. These different instances of unitary transforma- 
tions result from the irrelevance of the global phase in 
quantum mechanics and can be described mathemati- 
cally by multiplying with elements of the (finite) center 
of G. The center consists of those elements which com- 
mute with all elements of G. To find the time-optimal 
control algorithm, we may have to consider multiplying 
with different elements of the center. 

As exp(i7r/4)CNOT[2, 1] is an element of K, it takes 
negligible time to synthesize CNOT[2, 1]. In strong con- 
trast, exp(i7r/4)CNOT[l, 2] is not contained in K. Using 
the KAK decomposition, both exp(i7r/4)GNOT[l, 2] and 
exp(i7r/4)SWAP correspond to the same generator of A, 
given by n(—iS^ ' I x ) + 0(~iS a I x ), and the minimum time 
to synthesize each of them is equal to t m i„ = 7r. This is 
still the optimal time if we consider to multiply with dif- 
ferent elements of the center. 

We explicitly state the control algorithms: The unitary 
transformation exp(wr/4)CNOT[l, 2] is given by 

exp(i7rS' 2 /2) exp(«7r/ z ) exp(— iTrS a I x ) exp(— inl z ) 
= exp(inS,/2)exp(-il/H /J)&q? [-in H a (%) / 'wf] , 

where t' = —nJ/w^, mod 2tt > 0. Similarly, the unitary 
transformation exp(i7r/4)SWAP is given by 

e W 4 CNOT[2, i] e "/ 4 CNOT[l, 2] e -" r/4 CNOT[2, 1] 

_ e inSj2 e -inS a; /2 e -i3iTH /(2J) e iTrS y /2 e -it'H /J 



: exp [— iTrH a (ir)/wl] e 



I] p -iivS x /2 -iirH /(2J) 



TS y /2 



The corresponding pulse sequences are given in Fig. [3] 



V. PROOF OF TIME-OPTIMALITY 

In this section, we prove the time-optimality of the 
given control algorithms in order to synthesize unitary 
transformations in coupled fast and slow qubit system. 




FIG. 3: The figure shows the pulse sequences for synthesiz- 
ing the unitary transformations (a) exp(i7r/4)CNOT[l, 2] and 
(b) exp(i7r/4)SWAP, where R 5 = exp{iirS z /2) exp(-mS x /2). 
Since 1/J <C the length of the time intervals tt/w^ is 

larger as depicted. Refer to the text for details. 



A. The simple case 

All control algorithms, synthesizing a unitary transfor- 
mation in time t = Y]j tj , can be written in the form 

K' n+l expHC^WvOK ' ' ' K' 2 eM-it'iH (ipi)]K[, 

(15) 

where Kj G K take negligible time to be synthesized as 
compared to the evolution under H@ , tj , ipj €E K, and 
tj = tj/wj,. We can rewrite Eq. (fTS"]) as 

K n+1 ex V [-it n SPl x ]K n ...K 2 expl-ihS^Ku (16) 

where Kj G K. Equation (| L 6|) can be rewritten as 

K n+ i exp(p„) • ■ • exp(pi), (17) 

where pj — Kj{—itjS^I x )Kj and Kj are suitable ele- 
ments of K. Observe that the elements pj are contained 
in p. This follows from the Campbell-Baker-Hausdorff 
formula (see, e.g., Appendix B.4 of Ref. |44() and the 
fact that [t,p] G p (see Eg. [TO]). It was shown in Ref. [i| 
that for all time-optimal control algorithms the elements 
Kj can be chosen such that all pj commute. Therefore, 
all pj belong to a maximal Abelian subalgebra inside p, 
and we can find one K G K such that KopjKg 1 G a for 
all j. Using this result and results of Eq. ([50)) below, we 
can rewrite Eq. (I17p in the form 

K 2 ex.p(t n pn)---exp(tiPi)Ki, (18) 

where Pj = /? i (-iS ,/3 / 2; ) + a.j(-iS a I x ), 

ajf G {(-1, 0) T , (1, Of, (0, -1) T , (0, 1) T }, 



As expected, the maximal amplitude 
determines the optimal time. 



LJ,. 



(see Eq. ©) 



and Ki, K 2 G K. Equation p8|) can be simplified to 

K 2 exp[0(-iS I x ) + a{-iS a I x )]K 1 , (19) 

where a = J2j a j^j an d = YljPjtj- Assume that 
the unitary transformation to be synthesized is given 
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by one of its KAK decompositions K4 exp[ai(—iS^I x ) + 
a 2 (~iS a I x )]K 3 , where a,j € K and K 3 ,K A G K. We re- 
mark that the KAK decomposition is not unique, and 
we prove in Sec. IA 21 that different KAK decomposi- 
tions Kq exp[a[(— iS"l x ) + a' 2 (—iS a I x )]Kz l correspond 
to all values a'j = a,j + 2irZj, where Zj G Z and 
K$,Kq 6 K. We can choose ai and a 2 as those 
values of a[ and a' 2 such that |ai| + 1 0.2 1 is minimal. 
If I cii| + 1 02 1 > t, we cannot synthesize the unitary 
transformation in time t since all time-optimal con- 
trol algorithms are equal to Eq. (|19[) and \a\ + \j3\ = 

IEj + lE^I < E,(l«,-| + l&l)*, = Ei*i = *• 

For I ai| + |a2 < i, we can use the control algorithm 
exp(— iaiS^I x ) exp(— i-!rS x ) exp(— ia 2 S^I x ) exp(iirS x ) to 
synthesize the unitary transformation in time ai| + \a 2 \. 



B. The general case 

Until now, we have assumed that in Eq. tof, = 
u>i—J, i.e., we irradiate on the transition aa <-> a/3. More 
generally, we can irradiate on both transitions aa <-> a(3 
and (3a <-> /3/3. Hence, we substitute Eq. (fj]) by 

Hf h {t') = 2u I r (t'){b 2 cos[(wj - J)i' + 02(0] 
+ 61 cos[(wj + J)t' + 4> 1 {t')]}I x , 

where I&1I + I&2I < 1 (this ensures that the peak amplitude 
is 2ojI). Transforming into a double rotating frame by 

U lah (t') = exp(-it'uJsS z ) exp[-it'(wi + J2I z S z )]U Iot (t'), 

the evolution under the control Hamiltonian for time 
t' (with constant cuf, , $1 , <j>2 G K) generates a net ro- 
tation ZT('exp[— it'ul(bxSPl p + b 2 S a I q )}, where Z p = 
Z x cos(0i) + I y sin(0i), 7g = I x cos(0 2 ) + Zj, sin(0 2 ), and 
K'l G K. This can be rewritten as K[ exp(— itb)K' 2 , 
where b = bx(-iS I x ) + b 2 (-iS a I x ) G o, K[,K' 2 G K, 
and t = t'uif,. Therefore, any control algorithm generates 
in time t, a transformation (written as in Eq. (jTHJ) ) 

K n+1 exp[-it n b]K„ ■ ■■ K 2 exp[—itib]Ki, 

where tj is given in units of and Ei*i = ^ This 
generalizes the case of b\ = 1 and 62 = 0, treated in 
Sec. El] 

Similarly as in Sec. IV Al we obtain time-optimal con- 
trol algorithms as in Eq. (p~7j) . where pj — KjbK~ 1 and 

Kj are suitable elements of K. Therefore, Eq. (JT7J) can be 
transformed to Eq. I|18p . where the commuting elements 
Pj — fiji—i-SP I x ) + aj{—iS a I x ) are contained in the Weyl 
orbit W(b) = {KbK- 1 : K G K} n a, i.e., (/3 j> a j ) T is an 
element of the set 

{(&!, b 2 ) T , (bl, -b 2 f, (~h, b 2 ) T , (-h, -b 2 ) T , 

(b 2 , h) T , (-62, &i) T , (b 2 , -&i) T , (-62, ~6i) T } • (20) 



The Weyl orbit is induced by the map (K, b) 1— » KbK 1 , 
where & G a and the elements K G K are 

|id4, exp(— iirS a I z ), exp(— iirS^I z ), exp(— ml z ), 
exp(—iirS x ), exp(—iirS x ) exp(—iTrS a I z ), 
exp(—iirS x ) exp(—iTrS^I z ), exp(— i2irS x I z )^. (21) 

As before, Equation (|18|) can be simplified to Eq. (|19p . 
and we obtain \a\ + \(3\ < t(\bi\ + \b 2 \). Further- 
more, max{|a|, \/3\} < imax{|&i|, \b 2 \} holds. When the 
pairs (ai,a2) T and (bi,b 2 ) T satisfy max{|ai|, \a 2 \} < 
max{|&i|, I&2I} and |ai| + |a 2 | < |6i| + \b 2 \, then we 
say (ai,Q,2) T is r-majorized (6i,&2) T , i.e., (ai,a 2 ) T -< r 
(bi,b 2 ) T . The notion of r-majorization is equivalent to 
the condition that one element of a is contained in the 
convex closure of the Weyl orbit of another one (for a 
proof see Appendix I A ip . 

Given any unitary transformation G G G, let t opt be 
the smallest possible time such that 

(ai,a 2 ) T ^ r t opt (6i,6 2 ) T (22) 

and G = K 2 exp[ai(-iS^I x ) + a 2 {-iS a I x )]Ki with K 3 G 
K. Again, the KAK decomposition is not unique, and 
different KAK decompositions correspond to all values 
a'j = a,j + 2irzj, where Zj G Z (see Sec. IA 21 . Let 
us choose a,j as an element of [— tt, tt]. We prove in 
Appendix IA 31 that for such a choice of a,j, the equa- 
tion (&i,ci2) T -<r (o-i, a 2) T + 2ir(zi, z 2 ) T holds for all 
zi,z 2 G Z. This implies that the smallest t opt in Eq. (|2"2"1) 
can be achieved for ai,a 2 G [— tt,tt]. 

Then G cannot be synthesized in time t less than i pt> 
as for such a control algorithm the equation (a,/3) T -< r 
t(bi,b 2 ) T would hold, and this would contradict the 
minimality of t opt . In addition, G can be synthesized 
in time t greater than or equal t opt : It follows from 
(ai,&2) T -< r t op t(bi,b 2 ) T that (ai,a 2 ) T is contained in 
the convex closure of the Weyl orbit of i pt(fri,&2) T (see 
Appendix I A 1[) and we can synthesize G by convex com- 
binations of elements of the Weyl orbit of t op t(b\, b 2 ) T . 

Remark 3. Note, since (61,62) -< r (1, 0), it follows that 
the minimum time to produce any unitary transforma- 
tion can be obtained when all rf-amplitude is used to 
irradiate only on one nuclear transition (say aa <-> a(3 as 
in Fig. [T]) as described earlier. This justifies our initial 
choice of irradiating only on one nuclear transition. 

VI. CONCLUSION 

In this paper, we presented time-optimal control al- 
gorithms to synthesize arbitrary unitary transformations 
for coupled fast and slow qubit system. These control al- 
gorithms are applicable to electron-nuclear spin systems 
in pulsed EPR experiments at high fields. Explicit ex- 
amples were given for CNOT and SWAP. Recently, con- 
trollability results have appeared for c oup led electron- 
nuclear spin systems at low fields [13, where it is 
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shown that it is possible to synthesize any unitary trans- 
formation on the electron-spin system by only manipulat- 
ing the electron. New methods need to be developed to 
obtain time-optimal control algorithms in these settings. 
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APPENDIX A: PROOFS 
1. Convex closure of Weyl orbits 

Assume that a — a\{—iS^I x ) + a2(—iS a I x ) and b = 
bi(—iS^I x ) + b 2 (—iS a I x ) are elements of a. We prove 
that (ai,a 2 ) T is contained in the convex closure of the 
Weyl orbit of (h,b 2 ) T iff (ai,a 2 ) T < r (bi,b 2 ) T . 

Suppose (ai,a 2 ) T is contained in the convex closure 
of the Weyl orbit of (6i,6 2 ) T . Assume that |6i| > |6 2 |. 
Then, {a 1 ,a 2 ) T = TV w j( b j,U bj,2) T , where (6^1,^,2) 
belongs to the set in Eq. ([20")) (wj > and TV w j = !)■ 
It follows that l&^ij < \bi\ and \bj. 2 \ < \bi\. Therefore, 
I a,i I < |6i| and \a 2 \ < \b%\, implying max{|ai|, \a 2 \} < 
max{|&i|, \b 2 \}. Also note, |<zi| + \a 2 \ < J2j w j(\bj,i\ + 

M) = M + M- 

Suppose that (ai,a 2 ) T -< r (bi,b 2 ) T . The conditions 
max{|ai|, \a 2 \} < max{|6i|, |6 2 |} and |«x j + \a 2 \ < \h\ + 
\b 2 \ are equivalent to (|ai|,|a2|) T being weakly subma- 
jorized by (|&i|, |6 2 |) T - Thus, we obtain from Prop. 4.C.2. 
of Ref. [591 that (|ai|, |a 2 |) T = 

e 1 (|6 1 |,|6 2 |) T + e 2 (|6 2 |,|6 1 |) T + e 3 (|6i|,0) T 
+e 4 (0, \h\) T + e 5 (M,0) T + e 6 (0, |6 2 |) T = 
/i(N,N) T + M\h\, \h\) T + fsQhl ~\b 2 \) T 
+UHH \h\) T + f 5 (\h\, -\h\f + f 6 (-\h\, |6 2 |) T , 

where ej > 0, J2j e j — A = e i + ( e 3 + e 6)/2, f 2 = 
e 2 + (e 4 + e 5 )/2, and f k = e k /2 for k e {3,4,5,6}. In 
particular, we have that fj > (for all j) and J2j fj = ^- 



It follows that (ai,a 2 ) T = (ei|eti|, e 2 |a 2 |) T = 

/i(e3&i, tib 2 ) T + / 2 (e5^ 2 , £6^i) T + f3(^bi,e & b 2 ) T 
+fi{(-9b 2 , £io&i) T + fa{eub2, ei 2 &i) T + Uiexsbi, ei 4 6 2 ) T , 

for appropriate choices of ej G {1,-1}. We conclude the 
proof by consulting Eq. ([SIT]) . A Lie-theoretic proof can 
be obtained by following Thm. 2 of Ref. • 



2. KAK decomposition for elements of A 

We prove that the elements exp(a') G A equal to 
K\ exp(a)-ftT 2 are given by the elements {a'i,a' 2 ) T = 
(ai,a 2 ) T + 27r(zi,z 2 ) T , where K } E K, a' = a'^-iS 13 I x ) + 
a' 2 (-iS a I x ), a = ai {-iS p I x ) + a 2 (-iS a I x ), and Zj E Z. 

We can choose a' as a' = K(a + tyK^ 1 , where K is 
an element of Eq. (|2"Tj) and k € {q g a| exp(g) e K} 
(cp. Ref. [3, Lemma 2, and Ref. [17|, Prop. 4). Using 
the ansatz exp[a'{(—iS a I x ) + a^—iS^Ix)] — id4, where 
a'{, a 2 £ M, we obtain that a'{, a 2 £ {4irz : z G Z}. It is a 
consequence of Thm. 8.5, Chap. VII, of Ref. [43[ that {q £ 
a| exp(g) £ K} is equal to {qi(— iS a I x ) + q 2 {—iS a I x ) : 
51,32 G {2nz: z E Z}}. This completes the proof. 
We remark that exp[27T2i(— iS^I x ) + 2TTz 2 {-iS a I x )] = 
expfirtZii-iSPlz) + 2Trz 2 (-iS a I z )} for all z 3 E Z, where 

2irz 1 {-iS l3 I z ) + 2TTZ 2 (-iS a I z ) E t 



3. Proof of a majorization relation 

We prove that {a\,a 2 ) T ~< r (ai,a 2 ) T + 2n(zi,z 2 ) T 
holds for all z\, z 2 E Z, if we assume that a±, a 2 E [— 7T, it]. 

As the case z\ = z 2 = is trivial, we assume that 
\zi\ > or \z 2 \ > 0. We obtain that max{|ai+27rzi|, \a 2 + 
27rz 2 |} > 2tt — tt = 7r > max{ \a\ \ , |a 2 1 }, and the first 
condition in the definition of r-majorization is satisfied. 
The second condition I ai + 27rzi 1 + 1 o 2 + 27rz 2 1 > |eti| + |a 2 | 
follows from the fact that \a,j + 2irzj\ > \cij\ is always 
true. In particular, this is trivial for Zj = and it is a 
consequence of \aj + 2irzj\ > |(|27TZj| — |aj|)| > n > \cij\ 
in all other cases. The result follows. 
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